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In this Letter, we show that the classical 5*0(2, 1) symmetry of a harmonically trapped Fermi 
gas in two dimensions is broken by quantum effects. The anomalous correction to the symmetry 
algebra is given by a two-body operator that is well known as the contact. Taking into account this 
modification, we are able to derive the virial theorem for the system and a universal relation for the 
pressure of a homogeneous gas. The existence of an undamped breathing mode is associated with 
the classical symmetry. We provide an estimate for the anomalous frequency shift of this oscillation 
at zero temperature and compare the result with a recent experiment by [E. Vogt et ah, Phys. Rev. 
Lett. 108, 070404 (2012)]. Discrepancies are attributed to finite temperature effects. 
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The symmetry of a classical theory can be broken in 
the corresponding quantum theory. This is known as a 
quantum anomaly. Probably the best-known example of 
an anomaly is the conformal anomaly: While the classi- 
cal theory only contains dimensionless bare parameters, 
the quantization process forces us to introduce a regu- 
lator scale so that divergent terms are canceled. The 
bare quantities and the regulator conspire in such a way 
that the renormalized quantities are dimensionful. This 
startling phenomenon is called dimensional transmuta- 
tion. It is ubiquitous in quantum field theory, examples 
include QED, QCD, and the Gross-Neveu model. 

A quantum anomaly occurs for a two-dimensional (2D) 
gas of two fermion species as well. In this case the 
anomaly breaks classical scale invariance, and is there- 
fore referred to as a scale or conformal anomaly. The 
Hamiltonian of the system can be written as 



H 



drx 



A 



(1) 



We can read off the engineering dimension of the fields 
and the coupling Ij: The fermion fields '!/'cr(x) (with a =t 
, \,) have a dimension of momentum whereas the bare 
coupling constant A is dimensionless. Since we work with 
a local Hamiltonian, the quantum field theory contains 
divergences which must be removed by a renormalization 
of A. This is done in such a way that the leading-order 
term in the effective range expansion of the quantum- 
mechanical s-wave scattering phase shift is reproduced. 
Hence, the Hamiltonian ([T]) is referred to as the zero- 
range model. The renormalized dimensionful parameter 
is the bound state energy Ei,. With a hard momentum 
space cutoff A, is related to A and the regulator scale 
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Euler's constant. As the cutoff diverges, A — > go, the 
bare coupling A(A) —0 so that E^, is kept fixed. 

Experimentally, a 2D Fermi gas can be created by trap- 
ping an ultracold quantum gas in a geometry that tightly 
confines the gas in one direction [4H6] . Such a system can 
be described by the following effective Hamiltonian: 



H 



d'^x — 



(3) 



where uq is the azimuthal trapping frequency. Although 
the trapping potential explicitly breaks scale invariance, 
the system still exhibits a symmetry on the classical 
level [7j. The symmetry group is SO {2, 1), the Lorentz 
group in 2D. The generators of the corresponding algebra 
are constructed from the generators of scale and special 
conformal transformations, 



D 



(Px Ximji{x) and C 



£x'^i.li>^{x), (4) 



as well as the Hamiltonian of the noninteracting system, 
Eq. ([l|. ii = —i{ip^di'tp — diil)'^il))/2m is the current op- 
erator. Wc can express the Hamiltonian of the trapped 
system ([3| as a linear combination of two generators: 
iJosc = H + ujqC. D, C, and H obey the commutation 
relations [HI [9| 

[D, H] = 2iH, [D, C] = -2iC, and [H, C] = -iD. (5) 

The commutators in Eq. ([5| hold on the classical level, 
where they are understood in terms of Poisson brackets 
[A,BU ^ Jd'^if^jfr - Pffi). There are quan- 
turn corrections to those relations, and one of the main 
results of this Letter is to derive the anomalous correction 
to the commutator [D,Hosc], Eq- (13 1. The operators 

-uj^C±iujoD)/uJo, 



Li/2 = (L+±L_)/4, where L± ^ (H ^ 



Eb 



A2 



majpe^T^; 



and the Hamiltonian — Hosc/2ujq form the so(2, 1) al- 
(2) gebra [^1,^2] = L3, [^2,^3] = -Li, and [^3,^1] = ^2- 
This is an example of a spectrum generating symme- 
We set h = 1. a2D is the scattering length and is try: L_|_ and L_ are raising and lowering operators 
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with commutators [i7osc,^±] — i.2uJoL± and — 
—AHosc/^Q- If L± acts on an eigenstate of the Hamiho- 
nian i/osc with energy we obtain an eigenstate with 
energy E ± 2ujq. The excitations of the ground state 
were identified with breathing modes [7]. The work in 
Ref. [T has been generaUzed to the unitary Fermi gas in 
three dimensions for which the symmetry persists on the 
quantum level [SHTO]. 

Since the 5*0(2, 1) symmetry on the classical level is 
used to predict the breathing mode spectrum, the quan- 
tum anomaly should induce a shift in that spectrum. We 
provide an estimate for this shift in the hydrodynamic 
regime at zero temperature based on a recent Monte 
Carlo simulation of the equation of state [11. 

In this Letter, we show that quantum effects deform 
the algebra ([S]), which breaks the symmetry. We find 
that this anomaly is closely related to a two-body opera- 
tor known as the contact, which plays an important role 
in the description of various properties of an interacting 
Fermi gas. This work provides a new interpretation of 
the contact in 2D. The contact operator is defined as 

where d is the space dimension. As is well known, the 
matrix elements of / are finite [12] ■ They set the magni- 
tude of a class of exact relations that are referred to as 
universal or Tan relations [T3] . They include several ther- 
modynamic relations as well as the asymptotic form of 
various correlators, such as the momentum distribution, 
structure factors, or radio-frequency transition rates. For 
a review, see Ref. [Uj . Universal relations for Fermi gases 
in 2D have been derived in Refs. [TSHU]. Although we 
consider a fermionic theory in this Letter, all our results 
can be extended to bosonic systems. 

Quantum anomaly. — We start by deriving the form of 
the quantum anomaly. We show that the commutator of 
D and H receives an anomalous correction which is pro- 
portional to the contact operator /. This operator is not 
an element of the original so(2, 1) algebra: The algebra 
is deformed and the symmetry is broken. The effect of 
a quantum anomaly on a Bose gas at zero temperature 
has been studied previously in Refs. pOll^ . 

The change of the Hamiltonian under an infinitesimal 
scale transformation is given by [D^Hqsc]- We can sim- 
plify this expression using the Euler equation 

-im[ji(x), iJosc] = -djUij - x^TOWoViVvl^)' (7) 
where 11^ is the stress tensor of the zero-range model. 



n,,(x) 



2m 



(8) 



This yields, after integrating by parts. 



,{^)-2iLolC 



2iH - 2iu}, 



(9) 



If the last term on the right-hand side of Eq. ([9| were 
zero, the theory would obey the commutation relations 
listed in Eq. ([s]). In the following, we demonstrate that 
this term does not vanish and the trace of the stress ten- 
sor receives an anomalous correction. We choose a renor- 
malization scheme in which A^^ diverges as l/{d — 2) 
which cancels the prefactor in Eq. 

In order to determine the dependence of A on the reg- 
ulator, it is sufficient to calculate few-body matrix ele- 
ments. To this end, consider the scattering amplitude of 
two particles of opposite spin in the center of mass frame. 
The matrix element is obtained by summing a geometric 
series of "ladder" diagrams (see, for example, Ref. [TU]): 



A{E) = [m/\-GgMO)] 



(10) 



G^j.(r', r) is the bare propagator of a particle with energy 
E and reduced mass to/ 2. It is given by [10] 



gL(o,o) 



Q \47rsinhwo^ 



d/2 



_E d 

' 2wo 4 



1 



(11) 



Equation (111 diverges for d = 2. In order to extract 
the divergent part, we analytically continue d = 2 — e. 
To keep the coupling dimensionless, we introduce a di- 
mcnsionful scale in Eq. ([T]), which we identify with the 



inverse scattering length: A 



^2D 



A. This yields 



TO 

47r 



(0,0) 
2 

e 



In 



mL0Qa2Dt 
2^ 



-00 



E 
2ujQ 



(12) 



■00 is the digamma function. We choose a modified min- 
imal subtraction scheme with A~^(£) = {—2/e — 7b — 
ln7r)/47r. The poles of the scattering amplitude give the 
equation for the two-particle spectrum [22] . This gives 
rise to the relation 



[D, Hose] - 2iH - 2iu:lC 



2'Km 



(13) 



or, equivalently, \D,}i\ = 2iH + il /2Trm. The repeated 
action of H and D on / generates two-particle operators 
that contain additional derivatives. More generally, the 
correction term can be written as ijSxdH / d\, where f3\ = 
A^/27r is the beta function of the theory. The anomaly is 
due to the renormalization of the contact interaction and 
independent of the particular choice of the regulator. 
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Using an analytic result, the anomalous correction 
to a "hydrodynamic commutator" of a low-density 
Bose gas was determined in Ref. [20] . We reproduce 
this case by using the "adiabatic Tan relation" / = 
ATTma2BdH/da2B PS] in Eq. ^ (Note that (3x = X^/^tt 
for bosons). 

We can repeat the analysis for an untrapped system. 
Here, the system is not only 50(2,1) but also scale- 
invariant and the deformation of the algebra indicates 
a scale anomaly. Using the propagator of a homogeneous 
system, 0^0,0) = - j;^ dte(^+^°^^' {m/Antf^ in 



Eq. (|10|), we find: 
,(0,0) 



'^2D^homl 



m 
in 



2 , mEa^r.e^'^ 

-- +ln 

e 47r 



(14) 



The choice for A(e) remains unchanged. The trace of the 
stress tensor becomes 



1 



27rm 



(15) 



The conformal anomaly in the homogeneous case has 
been discussed in detail, both in quantum field theory [23] 
and as a quantum- mechanical toy model [U [M] US] . For 
the discussion of an enlarged model with an additional 
long-range interaction, see Ref. [5B]. 

Virial theorem and pressure relation. — Although the 
so(2, 1) algebra is deformed, we can use the new com- 



mutation relation ( 13 ) to derive one of the Tan relations, 
the virial theorem. The virial theorem relates the ground 
state energy Eq = (i?osc) of a many-particle system to 
the trapping energy w§(C). 

The thermal expectation value of 
an operator B is defined as {B) = 
tr{exp[-/3(77osc - M.iV,)]6}/tr{exp[-/3 {H^,, - fi^N,)]}, 
where fx^ is the chemical potential for the two spin 
species. Since [Hosc,N„] = 0, the thermal weight factor 
commutes with Hose and we have {[Hosc,B]) = 0. Setting 
B = D we can express the ground state energy in the 
trap as follows: 

Eo = {H + i^lC) = 2ul (C) - ^ (/) . (16) 

This form of the virial theorem has been obtained before 
in Refs. [151 [TH] using a different argument. 

A second important thermodynamic Tan relation is the 
pressure relation. It holds for a homogeneous system and 
links the pressure P and the energy density £. The pres- 
sure relation follows directly from the anomalous trace of 
the stress tensor, Eq. ( [Ts] ), which can be related to the 
pressure via 2PV = J aFx (nii(x)). Combining this with 
Eq. (15) gives the desired result: 



p = e 



(17) 



47rm ' 

where E and X denote the energy and the contact den- 
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FIG. 1; Anomalous frequency shift Suj at T = as determined 
from the fit to the equation of state in [11] . The red error band 
indicates the propagated error from the Monte Carlo calcula- 
tion. The plot also shows the exact results in the BCS (green 
dotted line) and BEC limit (red dashed line). For compari- 
son, we include experimental values for weak breathing mode 
excitations at T = 0.4:2Tp as reported in Ref. [30]. It re- 
mains an open question how thermal fluctuations affect the 
zero temperature results. 



as the virial theorem of a gas trapped in a box, but 
boundary terms have to be taken into account, so that 
{[D.H]) = 2iPV [37]. Again, the contact term in 
Eq. ( |17| ) is a direct consequence of the quantum anomaly. 
For a scale-invariant system, the pressure relation takes 
the well-known form P = E . 

Breathing mode. — As outlined in the Introduction, the 
50(2,1) symmetry implies the existence of hydrody- 
namic breathing mode excitations, the energy levels of 
which are spaced by 2u;o- A shift in the mode fre- 
quency is a manifestation of the quantum anomaly. A 
hydrodynamic description applies whenever the system 
is only slightly perturbed from thermodynamic equilib- 
rium which is maintained by frequent collisions of the 
atoms. In a strongly interacting regime where In fci?a2D ~ 
[3l|TT], the scattering cross section and, thus, the col- 
lision rate is large and we expect a hydrodynamic de- 
scription to be applicable. The effect of anomalies on 
the hydrodynamic equations has been much discussed in 
relativistic hydrodynamics, cf., for example, Ref. |28j . 

For a scale-invariant system, dimensional analysis dic- 
tates a quadratic dependence of the pressure on the den- 
sity n, P ^ . If we approximate the equation of state 
by a polytrope, P ~ n"'^^, where 7 is the polytropic in- 
dex, the linearized hydrodynamic equations can be solved 
(cf., for example, Ref. [29]). The frequency of the breath- 
ing mode in 2D is given by w^/wg = 27 -|- 2. We fit 
the numerical results for the energy per particle E/N in 
Ref. [TT] to determine the polytropic index 



7 ^ 



n dP 
P&n 



1 



e + g'+g74 



(18) 



sity, respectively. Equation (171 can also be derived where we parametrize i?/A^ = j^^(lnfci?a2D) with fc^? 



4 



'27m and define ^' 



94 



In Eq. 



18), the bound 
le result of this 



d In kfa2u 

state part of E/N does not contribute. T' 
calculation is shown in Fig. [l] The effect of the anomaly 
should be significant as we enter the strongly interacting 
regime, hikpa2Y) ~ 0. In this limit, we predict a shift of 
the order of 10% of the classical value. 

In a recent experiment by E. Vogt et al., the breath- 
ing mode frequency has been measured at much higher 
temperature, T/Tp w 0.4 30]. The authors do not ob- 
serve a significant shift from the classical value and ob- 
tain a result that is roughly consistent with the mean 
field prediction. This suggests that at finite temperature 
the effect of the anomaly is washed out by thermal fluc- 
tuations. It would be very interesting to determine the 
equation of state at flnite temperature to address this 
question further. (We checked that the theoretical re- 
sults in Fig. [T] are not significantly distorted for a slightly 
anisotropic trap and that they are robust under a vari- 
ation of the numerical data in [TJ.) The experimental 
results in Ref. [30] have been addressed in Refs. [3T1I32]. 
in which the damping of collective modes is analyzed us- 
ing kinetic theory. 

For large values of |lnA:Fa2D|, the energy per parti- 
cle can be stated in closed analytical form. We can use 
this to derive analytical expressions for the frequency 
shifts. In the limit lnA:i?a2D ^ 0, the system forms a 
Bardeen-Cooper-Schrieffer (BCS) supcrfluid and the fre- 
quency shift is 



tor and free Hamiltonian are 



1 



K 

2r] 



3 ' O , 



(19) 



where t] — lnkpa2D and k = 0.06 ± 0.02 has been de- 



termined in Ref. [TT]. Equation (19) is indicated by the 
green dotted line in Fig. [l] In the Bose-Einstein con- 
densate (BEC) limit, the system can be described as a 
gas of bosons with an effective dimer scattering length 
ttd ~ 0.55a2D [HI 133] • We obtain the anomalous fre- 
quency shift 



— = --+0(77-2 In 7;2), 
u}Q 4r/ 



(20) 



which is indicated by the red dashed line in Fig. [T] 
Equation (20) holds for a Bose gas as well, where Ud 



has to be replaced by the 2D scattering length 020- Fol- 
lowing the notation in Ref. |20j, we relate a2D to the 3D 
scattering length 030 via a^pn = ae^^^'^ /ne^'^'^^^ , where 
a = TTe^'"^+^{ C2D)^n &l, Czd ~ 1-47 and e = 030/^^0^, 
and = ^y2JrrlUJ^ is the harmonic oscillator length 
[5]. If we assume e <C min(l, 1/| Ino-j), the leading- 



order term in Eq. ( 20 ) reproduces the result for the 



anomalous frequency shift obtained by Olshanii et al., 
S = Suj/2ujo = a3u/4:V^az [20] . 

It is instructive to compare the present analysis to the 
Fermi gas in one and three dimensions (1D/3D). There, 
the commutator relations between the dilatation opera- 



ID 



[D, H] = 2iH 



«aiD 
2m 



I and 



3D : [D, H] = 2iH + 

47r rn 



(21) 
(22) 



As for the 2D case, we can use the commutator relations 
to derive the virial theorem and the pressure relation in 
ID [13 mi and in 3D [3 USUI]. Moreover, the sym- 
metry is broken explicitly at finite scattering length, and 
there is a parameter — the (inverse) scattering length 
fliD and aj^Q — that sets the strength of the breaking. 
For small parameter values, we can treat the correction 
as a small perturbation. This has been exploited to cal- 
culate the shift in hydrodynamic mode frequencies close 
to the unitary limit in 3D [35l |36] . 

In conclusion, we studied the effect of quantum fiuc- 
tuations on the symmetry properties of a 2D Fermi gas 
in a harmonic trap. We showed that there is a quan- 
tum anomaly, i.e., that the Pitaevskii-Rosch symmetry 
which exists on a classical level is not a symmetry of 
the quantum system. The anomaly manifests itself in a 
deformation of the Lie algebra associated with the sym- 
metry group. We derived that the anomalous operator 
appearing in the commutator relations is the contact op- 
erator and used this result to present a field-theoretical 
derivation of two thermodynamic Tan relations, the virial 
theorem and the pressure relation. We extracted the 
anomalous frequency shift of the breathing mode from 
Monte Carlo simulations at T = and compared the re- 
sult to recent measurements of the mode frequency at 
finite temperature. The findings of this Letter underline 
the subtle role that is played by the contact operator in 
the physics of interacting quantum gases. We can inter- 
pret it as an anomalous operator that is introduced by 
quantum fluctuations. 
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